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ABSTRACT 


A methodology for solar energy system design improvement has been 
developed and coupled with a constrained function optimization code 
resulting in an automated solar energy system design procedure. The 
scope of the methodology is limited to systems using flat plate 
collectors and water as the working fluid. 

Eight parameters have been included as independent design variables. 
The design variables included collector area, collector tile angle, 
collector and storage fluid stream velocities, and collector to storage 
heat exchanger dimensions. The procedure includes an accounting for 
economic parameters as an intimate part of the design process. The 
resulting methodology has been used for the design of solar energy 
systems which would use shelf item collectors for the purposes of 
determining the optimum design variable vector for a given situation. 
The methodology could also be used on a limited basis for collector desicn 
optimization by exploring the effects of changing selected collector 
parameters on system performance. The methodology is coded in the 
FORTRAN computer language under the name SOLOAD-1 (SOLAR ENERGY 
OPTIMIZATION ANALYSIS OR DESIGN) . 

Initial system trials indicate complete stability with minimal 
constraint activations. Based on the results of approximately fifty 
design experiments using SOLOQAD-1, new findings concerning optimm 
collector tilt angle and an invariant optimm collector flew factor 


are suggested. 
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NOMENCLATURE 


English Letter Symbols 


load heat transfer surface area, ft" 
collector area, et? 

Gaeaete,,  beu/ hr F 

collector cost, s/ete 

EUPLECose,. S/BtU 


systsm initial cost per unit collector area, 


o/ft 


e 


system installation cost per unit collector 
area, s/fite 


system lifecycle cost, $ 


system annual es power cost per unit 


collector area, o/fte yr 


Sycsteonetirecyele Cost per unit collector 
area, $/ft2 


2 
storage tank cost per unit collector area, S/Et a 


lifecycle operation and maintenance cost per 
unit collector area, o/ft* 


collector to storage heat exchanger cost per 
UM reCOLleCCtOr area, S/it* 


Specuelie heat, Bey/lb F 


monthly diffuse solar insolation on a horizontal 
surface, Btu/ft? mon 


collector loop inner diameter, ft 
collector loop out diameter, ft 
heat exchanger annulus inner diameter, ft 


value of the solar energy, $ 











D storage 


a) 


lifecycle cost, $ 
net Savings, $ 


annual fraction of the energy load provided 
by solar energy 


monthly fraction of the energy load provided 
by solar energy 


Coilecrormrecrtiiczrency Lactor 
collector flow factor 

system maintenance cost factor 
collector heat removal factor 
collector to heat exchanger efficiency factor 
present worth factor 
objective function 


eolteecror loop Flow rate per unit collector 
area, gpm/ft@ 


inequality constraint 


Soller insolationwon a Horizontal surface, 
Bem/nr EC 


extraterrestrial solar insolation on a 
hOEmZOnwaly Sumrace,, BCU/AL me 


solar insolation on a surface tilted towards 
the equator, Btu/hr ft 


Coneror tlag Lor COPES 
Neat exchanger length 

mass flow rate, lb/hr 
Humber CkecOnStmatnts 
number of design variables 
net present value, $ 


Mteration number 


rate of utilization of auxiliary energy, Btu/hr 





Mate nomena Loy thanster at the heat exchanger, 
Bew Ane 


rate of energy transfer from the load to the 
environment, Btu/hr 


GaeesOnmisetul energy collection, Btu/hr 
Sculceroructle COLELeECEIOn, Lactor 

PiguomOoremoncily GlEfuse ansolation on a tilted 
Sumemaces £o, Che monthiy direct insolation on a 


Homezental surtace 


average ratio of monthly direct insolation on 


a tilted surface to the monthly direct insolation 


on a horizontal surface 

ratio of monthly reflected radiation on a 
tilted surface to the monthly total radiation 
on a horizontal surface 

collector tilt angle, deg 

search direction 

ambient climatic temperature, F 

collector inlet temperature, F 

Soimececos OucLet EeMperature, F 

Storage tank outlet temperature, F 

eoulecceor loss Coefficient, Btu/hr oe F 
eoumlector loop tlow velocity, £t/hr 

Seepage Loop flow velocity, fEt/hr 

lower side constraint for i-th design variable 


upper side constraint for i-th design variable 


rate of work energy used to maintain system 
flow, Btu/hr 


vector of deisgn variables 
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Greek Letter Symbols 
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Q * 
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absorptance of the collector absorber surface 
move parameter in optimization problem 
declination angle of the earth, deg 

heat exchanger effectiveness 

lst collector flow parameter 

2nd collector flow parameter 

SQObllector efficiency 

unigue flow factor (optimum F") 

flow rate proportionality constant, gpm hr F/Btu 
geouna, GeElecclVity 

transmittance of collector covers 

latitude angle, deg 

Sunrise hour angle for horizontal surface, deg 


Sunrise hour angle for tilted surface, deg 
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ile RODUCT ION 


A. BACKGROUND 

The control and utilization of energy has become the 
major issue of this decade: the energy crises. The infla- 
tion rate of energy costs and particularly fossil fuel costs 
has the spector of uncontrollability as this decade comes to 
a close. Reference [1] reports a 30% increase in the price 
Of oil in 1978 alone. Reports from the media in early Septem- 
ber 1979 indicate the price of home heating oil has incresed 
70% since the 1978 heating season. The crisis center appears 
to be an emerging realization that the supply of available 
energy is in fact exhaustible. This realization has spread 
from a few to the masses; unquestionably the fruit of the 
OrganiZation of Petroleum Exporting Countries. 

In June of 1979, President Carter proposed to the Congress 
an energy strategy whereby solar energy would be providing 
20% of the Nation's energy by the year 2000. For the pur- 
poses of domestic hot water (DHW) and residential space 
heating this goal is technically achievable. The basic solar 
energy technical theory is well documented by Duffie and 
Beckman [2] and Kreith and Kreider [3] and the technology 
is continually being updated in the Solar Energy Journal [4]. 
The increased effort in solar energy research work is in 
clear evidence as the documentation in this journal appears 


to be on an exponential rise. 
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The treatment of solar energy system performance is 
typically a determination of the fraction of a given heat 
Meecdewhich is provided by solar energy. This fraction, f, 
is a nonlinear function of many system variables but typically 
increases in a monotonic fashion as a function of collector 
area. To simply seek maximum f and maximum collection of 
energy in an unconstrained fashion can result in extreme 
initial system costs. Accordingly, in order to assure a 
proper balance, technical performance and cost should be 
considered throughout the design process and economics should 
become an intimate part of the technical design problem; 
not an after thought at the end of system design. 

The economic ingredient in the problem of solar energy 
system design has received little attention in the literature. 
In view of this there is little information in the literature 
on how to proceed with deSign improvement of a solar energy 
system. Accordingly, a broad objective for this work was 
to develop a solar energy system model including economic 
considerations and to seek system design improvement by 
using the resulting eco~technical model together with an 


Optimization algorithm. 


Be SCOPE 

The conversion of solar energy into useful work covers 
a broad spectrum of collection schemes. These include thermal, 
photovoltaic, biological, wind, and ocean thermal energy 


conversion. Reference [5] 1S an excellent introductory 


aes 








information source covering many methods of solar energy 
conversion and has become a classic. This thesis effort 

was limited to the thermal conversion of solar radiation 

for the purpose of generating low temperature heat. Since 
low temperatures are involved, only flat plate collectors 
were studied. Flat plat collectors can provide temperatures 
on the order of 200° F compared to concentrating collectors 
which can provide temperatures as high as 6700° F. The 
restrictions on temperature were chosen in keeping with 
availability analysis; matching the solar conversion system 
to its task. The task of providing high entropy, low quality 
DHW and space heating energy 1S matched with low availability 
(1-€., minimum utilization of available energy in high 
quality, low entropy forms such as oil, gas, and central 
power). The collector working fluid was restricted to water. 
B@erder, since the collection of solar energy is a stochas- 
tic process, only analysis for long term performance was 
studied. The analysis of the dynamic performance of specific 
Systems in response to hourly climatic data has previously 
been conducted by Kline et al [6] and now forms the basis 


for long term analysis. 


See OBJECTIVE 

The objective of this thesis was to develop methodolo- 
eres for: 

1. System synthesis including a determination of optimum 


system design variables including collector area, collector 
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tilt angle, heat exchanger sizes, and flow rates. This 
analysis would be based on using shelf item collectors 
characterized by performance parameters Foto and FUL 
which have been determined by tests conducted in accordance 
with refs. [7 and 8]. 

2. Collector design optimization including a determin- 
ation of the optimum combination of collector design 
variables such as materials and geometries; and which 
would yield the optimum system performance. 

Fundamental to the effort was the utilization of: 

mec OPhs/CONMIN (Control Program for Engineering 
meeenesis with Constrained Function Minimization), a design 
improvement algorithm developed by Vanderplaats [9, 10, ll 
ama 12). 

2. Standardized long term solar energy load fraction 
correslations (f - CHART) developed by Kline [6]. 

3. Economic considerations combined with technical 
analysis. 

4. Long term climatic data obtained from (NOAA ) 


National Oceanic and Atmospheric Administration [13]. 


aS 








ooo new MODEL 


Peo VERALL SYSTEM DESCRIPTION 

A schematic diagram of the solar heating system is 
shown in Figure 1. The system consists of four heat 
exchangers, a storage tank and associated system piping. 
The first heat exchanger is called the collector and is 
used to transform incident solar radiation into useful 
thermal energy. At the second heat exchanger the collected 
energy is transferred to the secondary loop and stored in 
the form of sensible energy in the storage tank. The 
sensible energy 1s transferred from storage to the loads 
via the third and fourth heat exchangers. Energy transfer 
is achieved via four forced convection flow loops. This 
model was further simplified by assuming total utilization 
of the energy stored in the tank and thereby eliminating 


the need to detail the DHW and space heating heat exchangers. 


B. ENERGY BALANCE 


The steady state power balance for the system is 


+ > = 
QO W oe Q, (1) 
where, 

Q, = Useful rate of energy collection (Btu/hr) 

ae = Work energy utilized to maintain system flow 
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Q [eh acme Ue! i zaelOn OL auxiliary energy 


5 which is required in addition to solar 
energy in order to meet the DHW and space 
heating load 

Q, = Rate of energy loss by the system 


For the purposes of this thesis, it has been assumed 
that aes is a known or estimated parameter. Most analyses 
Mmaeene literature do not account for he Since it is typically 
assumed to be a very small fraction of the energy required. 
The percentage of the energy load which is supplied 


by solar energy is defined as: 


thi 
r 


where the summation is taken over 12 monthly values. 

Emisea monotonically increasing function of collector 
area. No optimum f is evident. Figure 2 depicts a typical 
f distribution resulting from a design problem. Since 
initial costs are strongly dependent on collector area (See 
Section III), attempts to seek maximum £ can result in very 
Peagieinitial system costs. At this point the motivation to 
include economic considerations as a part of the system 


design model become clear. 


Ce. ECONOMIC ANALYSIS 
The matter of economics has been treated as 1t evolved 
in this study: a simple common sense approach. This 


approach includes: 


Ly, 











1. An economic equation which utilizes technical 
parameters obtained from mass, momentum and energy balance 
considerations together with economic parameters. 

2. Standardization of the economic equation by char- 
acterization of all terms in present dollar values. 

The economic analysis 1S similar to an energy balance. 
Considering the economic system as the monetary account 
of a consumer and the time period as the economic lifetime 


of the energy conversion equipment: 


= - 
Din Pout D storage (3) 
where, 

D., = The value if dollars of the energy 
produced by the equipment 

D = The expense dollars to obtain, operate 

Sut 
and maintain the equipment 

D = The net savings 

storage 


The rational consumer seeks maximum D This 


storage’ 
methodology iS a Simple capital budgeting technique and is 


contained in any standard reference in accounting or finance. 


Sreeecie OBJECTIVE FUNCTION 

Implicit in the use of the economic model 1S an assump- 
tion that the terms in the equation can be expressed as 
Variables of the technical functions which have been 


defined by the physics of the engineering problem. For 
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the solar energy system utilized for domestic hot water 


and space heating the transformation is simple: 


Net Savings = Fuel Savings - System Costs 


This equation becomes the objective function by 


standardizing all dollar amounts to present values. The 


resulting equation is: 


NPY = £0,C Fi - ac. (4) 
where, 
NPV = Net present value of the solar investment, $ 
Ce = Fuel cost for the energy replaced by solar 


energy, $/Btu 


Bo = Present worth factor which standardizes fuel 
Savings during system life into present 
dollars, years 


A = Collector area, fee 


€ Seo VvSeeimiiteevalesGOSt per unit collector 
area, $/£t2 
The effort now proceeds to determining each of the 
above dependent parameters and to obtaining those combina- 
tions of these variables which will result in maximum 


values of NPV. 
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Pio cole OF FOUATIONS 


A. GENERAL 

The objective function has been fully developed by 
examining each of its parameters individually. Completion 
of this development has resulted in anon linear function 
of eight independent design variables; a hypersurface in 
eight space. When formulated in this manner, the problem 
fomelearly not amenable to analytical solution. Ina very 
Simple manner however the problem can be interfaced with 
the design improvement algorithm COPES/CONMIN (see Section IV). 
The primary contribution for this effort comes from Kline [6] 
who developed the f-chart correlations for determining monthly 
solar energy load fractions fi: The annual load fraction 
£ is then computed by a weighted average of the monthly 


values. 


pees BC LE TC” PARAMETERS 
1. Annual Energy.Load, QF 

The annual energy load consists of two basic ingredi- 
ents; the DHW load and the space heating load. Assuming an 
average mean ground temperature the month to month DHW load 
was considered constant except for the variation due to month 
length. The methods for determining the space heating loads 
typically follow the ASHRAE manual, reference [14]. Refer- 


ence [15] is recommended as a text for space heating load 
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computations. For the purposes of this study representative 
building loss parameters were used. Typical building con- 
ductance or (UA) | wares Ornso,000,, 20,000 and 10,000 
(Btu/deg F day) have been used based on a standard building 
@e 1/750 ee of floor area, a heat transfer surface of 5000 
ee and building loss coefficients, OS ee, .1/ 5 and .09 
pen /ft- hr F respectively. The conductance, UA, is the 
Space heating load at design conditions which has been 
estimated in the manner of Ref. [14] and divided by the 
design temperature difference. 
PeeeAnoual Solar Energy Load Fraction, f 

f is a function of collector performance parameters, 
thermal physical properties of the working fluids, fluid 
flow rates, heat exchanger performance parameters, collector 
merce angle, climatic conditions and latitude. The climatic 
conditions determine the load distribution (heating degree 
days), the solar energy flux distribution (insolation) and 
the environmental stress on the collector (ambient tempera- 
ture). 

a. Collector Performance 

A determination of collector performance proceeds 


From the well known collector equation of Hottel and Whillier 


[16] which is the result of an energy balance on the collector: 


Q = BU 


7 A_(IRta - UL (tT, ~ To)? (5) 


i 


Za 





where, 


One = Useful rate of energy collection (Btu/hr) 

P = Collector heat removal factor (dimensionless) 

UL -meclleccor Hoss coefficient (Btu/hr £t* F) 

AQ = Collector area (£t2) 

I = Solar insolation on a horizontal surface (Btu/hr £t*) 

R = Collector tilt correction factor (dimensionless) 

T = Transmittance of collector cover system (dimensionless) 

a = Absorptivity of the collector absorber plate 
(dimensionless) 

T. Secoklector fluid inlet temperature (F) 


Ta = Climatic ambient temperature (F) 


The analytical expressions for P and U, are complex 


L 


and result from a lengthy development which can be found in 
references [2 and 3]. The computation of nee and UL using 
these expressions is not necessary for system design analysis 
uSing shelf item collectors since the parameters FU, and 
FUta can be obtained from collector efficiency tests. 


Collector efficiency is defined as: 


QA 
IR 


Assuming that U. is not sensitive to changes 


L 
in the environment and combining equations (4) and (5) results 


in: 


inet= P to - PU, (Tt, - T.) / IR) (7) 


we 





Equation (7) is linear in the collector 
parameter (T. = T,)/7IR and forms the basis for determining 
collector performance. It is clear that the parameter 
PUL is obtained from the slope and Pto is obtained from 
the intercept of equation (6). Some typical curves of 
SPembector efficiency are included in Figure 3. 

Two standards have been developed for collector 
testing to measure collector performance, namely references 
mi) and [8]. The standards differ in the independent 
variables used in the performance characterization. 
Reference [7] the ASHRAE standard uses T = T. = Ti 
ref [8] the NBS standard uses T = (T + TS) 72. where Ty 
mememe collector outlet temperature. All collector 
performance test data used in this effort is based on the 
ASHRAE standard. All collector data used in this effort 
is based on actual collector tests conducted by NAVFAC 
and reported via ref [17]. 

Hee Collector to Storage Coupling 

The collector becomes coupled to the storage 
tank via the collector to storage heat exchanger following 
the method of de Winter [20]. It is assumed that the rate 


Of energy transfer at the heat exchanger is equal to the 


rate of useful energy collection: 
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where, 


T 
e) 


lk 
Ss 


Maximum system temperature or collector outlet 


Minimum system temperature or storage outlet. 


Equations (5) and (8) are combined to eliminate 


T. and T dependence to give: 


Q1 = FU, (IRta = U, tT. = T,)) (9) 
where, 
F. = ee (10) 
: -l 
EF U_A (me) 

| = {1 ee ee | (a) 

(mee) Sa eile 3 ee 

sian =: p min 


eee cOlar Insolation at Optimum Tilt Angle, T, 
The solar radiation intensity on a tilted 


Smetace is by definition: 


= JB dey) 


Solar insolation is typically measured and 
reported for a horizontal surface. Improvements in a 
Sema cCOliector installation are sought by tilting the 


collector to the optimum angle. The computation for 








optimum angle is iterative for each design situation. 
Several rules of thumb have evolved: 

(1) Collectors should be oriented at a slope 
of .9 times the latitude angle for maximum annual collection. 

(2) For DHW systems where loads vary little 
during the year the best angle of tilt 1S equal to the 
ieiitude angle. 

(3) The optimum collector tilt angle for least 
cost per Btu delivered for building heating is approximately 
the latitude angle plus 15 degrees. 

HiemMiecenOdeolLegy stem determining collector tilt 
angle used in this effort was as follows: 

(1) Horizontal monthly data was obtained for 
97 locations from NOAA in reference [13]. 

(2) An algorithm was developed to compute 
Memtnly horizontal extraterrestrial radiation intensities. 
ies Computation is a function of latitude, daily hour 
angles on horizontal and tilted surfaces, daily declination 
angle, and collector tilt angle. The azimuth angle was 
always chosen for a due south collector orientation. 

(3) The methodology follows that of Liu and 


Jordan in reference [19]: 


56 = ae 


HA. 


d 
T Ra + R (13) 


a ae 


D 


Pas) 





Z S 


C= 1.3903 - 4.0273() + 5.5315(4)° - 3.108(=)" (14) 
o) o) O 

w ‘sin(o-s)sin §6 + cos(g-s)cos § sin w_' 

d w Sin d Sin 6 + cos ¢ cos § Sin We 

We = minlcos ~(-tan(¢-s) tan one eoene (ten @ tan 6)] (iG) 

= oe * fiesta @? tan 6] (7) 

R, = + (1+ cos s) (18) 

d 2 
Ik 

R = 5 ee (COS Srp (19) 

where, 

R = monthly slope correction factor 

Rp = average ratio of monthly direct radiation on a 
trlted surface to monthly direct radiation on a 
horizontal surface 

Ra = ratio of monthly diffuse radiation on a tilted 
surface to monthly diffuse radiation on a 
horizontal surface 

fee ratio Of monthly reflected radiation on a tilted 

P Surface to monthly total radiation on a horizontal 

surface 

$ feeraclenot monthly diffuse radiation on a horizontal 
surface to monthly total radiation on a horizontal 
surface 

= Rotel Osoremonciliy total radiation on a horizontal 

oO Surface to monthly total extraterrestrial radiation 

on a horizontal surface 

w= Sunrise hour angle on a horizontal surface 

a Sunrise hour angle on a tilted surface 
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pee— collector tilt angle 


d = latitude angle at collector site 
BS = declination angle of the earth 
0 = ground reflectance for area adjacent to collector. 


3. Fuel Cost Parameters 
Fuel cost parameters included a selection from among 
Mme types (1.e. oil, electricity, or gas), unit of issue 
cost, fuel heating value and efficiency of the auxiliary 
heating system. The SOLOAD-1 system allows complete freedom 
in the selection of these parameters. A typical set of fuel 


cost parameters used in a deSign experiment included: 


Fuel: Qari: 

Gmat Cost: =e ($/Gal) 

Heating value: 427000 (Btu/Gal) 
Bit LCLency : Hey 

Resulting Fuel Cost: C. = $9.05 per ioe Bou 


Ys 


System Cost Parameters 





System cost considerations typically include 
initial costs and life cycle costs for operation and 
maintenance. Detailed guidance for complete cost consi- 
derations is included in ref. [3]. Further, reference [20] 
contains guidance to be used in determining costs for US 
Navy installations. The following cost parameters were 


used in SOLOAD-1: 


fan 





G = AC (20) 


where, 
CL —welOtal system life cycle cost 
Ce = Total cost per square foot 
E — er re (2k) 
Ss 5 V 
= + + + Ze 
C5 C Cc, +C.,+ ¢, (As 
E =e emp 16.) (23) 
V A eh 
where, 
C. = Initial system cost 
C., = Operation and maintenance costs 
Co = Collector cost 
Cy = System installation cost 
oon Storage tank costs 
oe = Collector to storage heat exchanger costs 
os = Pumping power costs 
a = Maintenance cost factor 
Py = Present worth factor 
., oN 
phos tc oe (24) 
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where, 


1 = i- 
cs Tass (25) 
i = Discount rate 
j = Fuel inflation rate 
N = System life in years 


It should be noted that Fe ©CGUmMswexpl@ertliy in 
mae fLairst term of the objective function and implicitly 
in the second term. It should also be noted that initial 


costs are not amortized in the SOLOAD-1l model. 
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IV. NUMERICAL OPTIMIZATION 


A. SIMPLE OPTIMIZATION 

Design problems typically seek the minimization or 
maximization Of an appropriate parameter within a frame- 
work of constraint specifications. The parameter to be 
optimized may be a function of several deSign variables and 
is termed the objective function. Other parameters which 
may be separate functions of deSign variables must not 
exceed specified bounds for the design to be acceptable. 
These parameters are termed design constraints and are not 
to be confused with limits which may be set on design 
variables which are usually termed side constraints. 

Engineering problems can be numerically coded for an 
analysis (once through) solution. The simplest scheme for 
Optimization may consist of a series of loops through the 
computer code which may cycle many combinations of design 
variables. The combination of variables which provides the 
best design and which also satisfies the ree is 
then considered optimum. This approach may be economical 
for a design problem with just a few design variables and 
short computer time requirements. A design problem with 3 
design variables, ten values for each design variable, and 
One-tenth seconds central processing unit (CPU) time for 
each analysis would take a total of 100 seconds of CPU time. 


The solar energy optimization design problem as characterized 
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by the code developed by this thesis has a minimum of 8 
design variables; each analysis of its objective function 
requires about 2 CPU seconds. Using the simple approach 
and assuming ten values for each design variable would 
result in a CPU time of 68 years for each design problem. 
Clearly, a more rational approach to optimization is 
necessary. 

Vanderplaats [9] suggested that many special algorithms 
for numerical optimization have been proposed in recent 
years, but that in many cases unsuspecting practitioners 
mime their particular optimization problem unsolved only 
after large amounts of time and effort are expended. This 
can occur usually because of inexperience by the practitioner 
in determining the limitations of specified algorithms. 
Vanderplaats [10, 1l, and 12] has developed a FORTRAN 
coded algorithm capable of optimizing a very wide class of 
eagancering problems. The system includes COPES (Control 
macgram fOr Engineering Synthesis) and CONMIN (Constrained 
momecrtOn Minimization). This optimization system is referred 


eeeas COPES/CONMIN. 


Poa COPES/CONMIN SYSTEM 


i Lermino logy 
CONMIN is a FORTRAN program in subroutine form for 


the solution of linear or non-linear constrainted optimiza- 


tion problems. The user prepares an analysis program. The 


ee 








program must be named SUBROUTINE ANALIZ. The process of 
computer aided design or of trade off studies with a minimum 
of man-machine interaction becomes fully automated via the 
COPES program. Three basic definitions are required: 

Design Variables. Those parameters which the 
Optimization program (CONMIN) is allowed to change in order 
to improve the design. Design variables appear only on the 
right hand side of equations in the analysis program (ANALIZ). 
Limits imposed on design variables are termed side constraints. 

Objective Function. Usually the single parameter 
which is to be minimized or maximized during optimization. 
The objective function always occurs on the left side of 
the equation in the analysis program. (Refer to [12] for 
exceptions.) The equation defining the objective function 
Maye be Linear or non-linear, implicit or explicit, but must 
be a function of the design variables to be meaningful. 

Design Constraints. Any parameter which must not 
exceed specified bounds for the design to be acceptable. 
Constraint parameters always appear on the left side of the 
Penetraint function equations. Constraint functions may be 
mimear OF non-linear, implicit or explicit, but must be 
functions of the design variables. 

Assuming that the optimization process requires the 
maximization of a particular objective function, the general 


Optimization problem can be stated as: 


BV 





a. Find: X which maximizes f (X) 
5 Besibleyj iayere “cee 
(1) CONSTRAINT EQUATIONS G, (X) 0 ge) NCON 
(2) SlwE CONSTRAINTS VLB 5 Ss Xs S VUB ., i 1 NV 
Where, X = X(Xy, Kos a, X,) is the vector 
Of NDV design variables, F(X) is the objective function and 
G, (x) are the set of NCON constraints. VLB . and VUB. are 
the lower and upper bounds respectively on each design 
variable. 
2. Methodology 

The solution process proceeds as follows: 

a. The user prepares an analysis subroutine which 
defines X, F(X) and G(X). This subroutine must be named 
ANALIZ. ANALIZ must have three segments; input, analysis 
m@emoucout keyed to COPES flags, ICALC = l, 2, or 3 
respectively. 

ce Lhe Wser prepares an input data file for COPES 
which includes a wide variety of system options, appro- 
priate matching mechanisms between ANALIZ and CONMIN and 
the constraint boundaries. 

ee. Using the initial vector of design variables 
COPES obtains an initial solution from ANALIZ and subsequent 
solutions by updating X as determined by CONMIN. Any analy- 
Sis solution which satisfies the constraint equations and 
the side constraints is a feasible design. If an analyti- 


cal solution violates any of these constraints it is an 


infeasible design. The minimum feasible design is optimal. 


Be 





The feasibility determination includes: 

(1) If a constraint equation is violated 
mere., if G(X) > 0) then the jth constraint is violated. 

(2) If a constraint equation equality 
condition is met (i.e., if G(X) = 0) then the jth con- 
Pemaimnt is active. 

(Sere oecCOnstraint equation condition 
Moremet (i1.e., if G(X) “Oe tense leweOnStEralnt 1s inactive. 

Note that CONMIN is designed to minimize objec- 
faye fLUunctions; the process of maximizing an objective 
function is concerned with minimizing the negative of an 
objective function. 

(4) Similarly, side constraints may be 
inactive or active but side constraints will never be 
violated in a particular analysis computation because 
they are specified limits not dependent variables which 
is the case for design constraints. 

(5) All inequality conditions are repre- 
Peieecwoy 2 band around the zerc condition due to computer 
bimitations in defining zero. 

d. If the initial analysis solution is infeasible 
CONMIN moves towards a feasible solution by adjusting the 
design variables appropriately. The optimization process 
then proceeds in an iterative fashion as follows: 

(1) The iterations are governed by the 


meemrsion relation: 
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where, 
q = iteration number 
a* = a scalar move parameter which defines the 
distance of travel in the direction of 
Search 
§ = direction of search 
(2) S$ is determined such that F(X) will 
be minimized without violating any constraints. CONMIN 


Smiemlates the gradient of F(X) by using finite difference 
techniques. Because no constraints are violated, the 
greatest improvement in F(X) will be realized by moving 

in the direction of steepest ascent, the gradient of F(X). 

Gye Once the warective is known, the 
move parameter, a*, which will allow the largest magnitude 
improvement in F(X) 1s to be found. A one dimensioned 
Beareh of the S direction is carried out until the best 
improvement number, a*, is found. 

(4) CONMIN utilizes methods other than 
the method of steepest descent in determining S particularly 
in the presence of active constraints. These methods 
include the method of conjugate directions developed by 
Fletcher and Reeves [21] and the method of feasible direc- 


tions developed by Zoutendijh [22] and implemented by 


Vanderplaats and Moses [23]. 
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e. CONMIN continues to iterate for an optimal 
design by computing successive S and a* values always 
keeping within the defined constraints. If there is no 
relative or absolute change in F(X) for three successive 


iterations, the optimum is considered found. 
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Veer DO LOAD] sy STEM 


A FORTRAN coded algorithm for the analysis of systems 
containing shelf item collectors has been developed. The 
algorithm has been named SOLAR ENERGY OPTIMIZATION ANALYSIS 
Seer oltGN (SOLOAD-1). The portions of the algorithm 
reserved for design optimization of collectors were not 
completed and have not been included in SOLOAD-1l. The 
interface between SOLOAD-1 and COPES/CONMIN has been 
developed via the Control Program (CP-67) and the Cambridge 
Monitor System (CMS) a time-sharing system developed for 
the IBM System 360 Model 67. SOLOAD-1 includes: 

A. Three executive programs to provide the interfaces 
among the system elements and to initiate design problems. 

B. An analysis subroutine which includes the system 
of equations covered in Section III. A summary of the 
objective function, design variables, and design constraints 
is included in Appendix A. 

C. Three auxiliary subroutines including a calendar 
data array, climatic data and user defined input parameters. 
D. Data files for user input to SOLOAD-1 and COPES/ 

CONMIN. 

The interface among SOLOAD-1 system elements and 
COPES/CONMIN is shown in Figure 3. SOLOAD-1 programs and 
Subroutines are included in Appendix B. Data files are 


included in Appendix C. 
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Vila RE Suis 


A. GENERAL 

The first objective of developing a methodology for 
the analysis of systems with previously defined collector 
design was achieved. The second cbjective of developing 
a methodology for optimizing the collector component design 
was not achieved. However, it appears that the potential 
exists for using the system analysis methodology to obtain 
improvements in collector design. For example, if cost 
data on various types of collector cover plates or collector 
absorber surfaces can be correlated with collector per- 
formance parameters, then the system model could be used 


for quick checks on collector design improvements. 


ee VERIFICATION 

A single design problem was selected and started using 
10 different sets of the starting design variable vector. 
The same optimal design result within + 0.5% was achieved 
for each run. This was the only actual verification deemed 
possible at this time due to lack of any known optimization 


data in the literature. 


fee eAPERIMENTS 
Approximately 50 design experiments were analyzed. 
Each problem was characterized by a unique identification 


number which could tie together the location, economic 
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environment, the collector, and the space heat transfer 
coefficient, UA, for the experiment. For each experiment 
an input parameter summary report and output parameter and 
results summary report was generated. These reports are 


included in Appendix D. 


Peo PECIFIC RESULTS 
ieee Constraints 
The only active design constraint in most experi- 
ments was the tube thickness in the collector loop at the 
heat exchanger. Some experiments resulted in no active 
design constraints. There were no active side constraints 
in any experiment. 
meee COLLlector Flow Factor 
The collector flow factor F" is defined as the 
ratio of the collector heat removal factor to the collector 


Peeirerency factor: 


The flow factor result obtained in all optimization experi- 
ments was the same. F" = .948 + .008. 
pee collector Tilt Angle 
Each experiment conducted for a particular geo- 
graphical location resulted in the same collector tilt 
angle. However, experiments for different locations but 
the same latitude resulted in significantly different 


Speimum collector tilt angles. 
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Vise, CONGZUSTIONS 


A. The SOLOAD-1 system in conjunction with COPES/CONMIN 
appears to offer the potential for further improvement and 
potentially a valuable automated technique for solar energy 
system design. 

Seeeetae rule of thumb” typically used for collector tilt 
Prgeree Optimization (1.e., latitude plus 15°) should be 
used with caution since preliminary results indicate a 
strong tendency for climatic dependency. 

C. The uniqueness of the flow factor was suggested by 

the continued result of .948 + .008 for each experiment. 
Pending confirmation by further testing it appears that 

a simple correlation for determining optimum collector 
loop flow rates may be available. Based on a flow factor 
of 0.948 the resulting correlation would be: 


G = .01955 PU (gpm/£t* Gtrecoltector abea) 


L 


This correlation is developed in Appendix E. 
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VIII. RECOMMENDATIONS 


A. The basic model of SOLOAD-1 should be upgraded to 
include specific characterization instead of simple 
parameter selection and input for: 


1. Pumping power in all four loops 


2. Inclusion of the DHW heat exchanger 
3. Inclusion of the space heat exchanger 
meeemnclusion Of building loss coefficients as a 


design variable. 
B. Prior to any additional experiments with the present 
model, a complete survey of the industry should be con- 
ducted for collector parameters and costs. Following this 
Survey, a series of experiments should be conducted in 
search of correaltions among collector performance parameters 
ee. 2., F(ta) and PUL) collector unit costs, and system 
performance. 
C. A larger sample size should be used to verify the 
uniqueness of the optimum collector flow factor (F") as 
suggested by the results of 50 experiments. 
D. The remainder of the NOAA climatic data bank (i1.e., 67 
more cities) should be included in SOLOAD-1l. The optimum 
eolltector tilt angle could then: be computed for the 97 
NOAA regions. 
E. The model should be upgraded to accommodate analysis 


for systems with air as the working fluid. 
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F. The model should be upgraded to accommodate amorti- 
zation of the installation cost instead of just initial 


cash payment. 
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meojke 3 Typical Collector Efficiency Curves 
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FIGURE 4 SOLOAD - COPES/CONMIN Interface 
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APPENDIX A 


SUBROUTINE ANALIZ SUMMARY 
Pee OBJECTIVE FUNCTION 
= £ ' = 
NPV f Qe C. Pe Aa ae 


B. DESIGN VARIABLES 


Collector Area A. 
Golleccorn, frie Angle S 
ColMwecerem Eoco, 1[.D. d. 
Colleceor Loop 0.D. 7 
Heat Exchanger I.D. ds 


Collector Flow Velocity Me 
Storage Flow Velocity V 


Heat Exchanger Length L 


=e DESIGN CONSTRAINTS 


G) = qe = 

oD 7 ora Ss q, 

G, = Reynolds number, Collector loop 

Gy = Reynolds number, Storage loop 

Ge = Capacity ratio, Ca oR 

Ge = lst Flow parameter, C1 = where oe UL 

Go = 2nd Flow parameter, C5 = Cee 
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| APPENDIX B 


SOLOAD-1l Computer Program 
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APPENDIX D 


EXPERIMENT REPORT SUMMARIES 


Each optimization problem or experiment is represented 
by an input summary report and an output summary report. 
These reports are tied together by a unique identification 


as follows: 


where, 


LL = location identification number in 
accordance with Appendix C 


Ny = present worth factor identification number 
iaees = £8.22 
1 
2- F! = 24.34 
a 
No -mce WLeCrEOnm TGCNEI: leation. number 


1 - Solarnetics 

2 - American Sun 

3 ~- Federal Prison Institute Double Glaze 
N = heat load conductance 


1 - UA 


30000 Btu/HDD 


2 = UA 20600 "Bru7 wbL 


Ww 
i 
CG 
ao 
il 


10000 Btu/HDD 
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APPENDIX E 


POTENTIAL CORRELATION FOR OPTIMUM COLLECTOR FLOW RATE 


Reference [2] defines the collector heat removal factor, 


F , as: 
r 





a T {l - exp ( Ge Pai (Bl) 


meee hner, Ehesecomlecceorettonw Eactor, F"”, is defined as, 





i Ge ~P'U, 
FE = Pr = ro lee CxXD ( Ge. es: (B2) 


Reference [2] shows that F" approaches the unity value 
asymptotically as the parameter CU ae increases. 
Reference [3] recommends capacity rates of 10-15 lb/hr ete 
meo2z23 ~ .0334 Gpm/ft*) as the best compromise among 
collector heat transfer coefficient, fluid pressure drop, 
and energy delivery. The unit ft. refers to collector area. 

The results of this thesis suggested a correlation 
between the collector performance parameter PU, and the 
Capacity rate. Once these results are verified by further 
testing including model changes to include fluid pressure 


drop, a simple correlation may be available and follows 


meom : 


228 





ifemecollows from 


Gr, 
* 

Bet; oy = 

Therefore, 

Or 

Or 





c<* = The optimum flow factor determined from 
computer experiments. 


F" = 


= Gc /F'U 
o/ 


= Be Bt 


t- 





v 
NO 


L 


L 


equations (B11) 


and (B2) 


sl exp (-1/z,) 


=in (l= 


that: 


The optimum parameter which corresponds to 


<* and is obtained by solving equation 


‘clepe oy 


| 


AN IN 
+1 4 
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(B2) 








Ot, 


: o4 
G = mF PUL /e 
Or, 
= KKP UJ 
G K oe L 


Based on the results of a limited number of computer 
experiments all yielding an apparent invarient flow factor, 
fom 948: 

aoe — 701955 (Gom hr F/Btu) 
Or, 


_ 2 
G = -O1955F U (Gpm/ft -) 


L 
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